#ifndef BERNSTEINBASIS_HPP
#define BERNSTEINBASIS_HPP

class BernsteinBasis
{
  public:
    BernsteinBasis( int p, double xi);
    ~BernsteinBasis();

    /* deri i=0,1,2 represents the function value,
    1st order derivative, and 2nd order derivative */
    double* der0;
    double* der1;
    double* der2;
    
    int getorder() const
    {return order;}

    double getIntPoint() const
    {return x;}
  private:
    void initialize();
    int order;
    double x;

};


#endif
